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Abstract. We study semi-classical limits of cigcnfunctions of a 
quantized linear hyperbolic automorphism of the torus ( "cat map" ) . 
For some values of Planck's constant, the spectrum of the quan- 
tized map has large degeneracies. Our first goal in this paper is 
to show that these degeneracies are coupled to the existence of 
quantum symmetries. There is a commutative group of unitary 
operators on the state-space which commute with the quantized 
map and therefore act on its eigenspaces. We call these "Hecke 
operators" , in analogy with the setting of the modular surface. 

We call the eigenstates of both the quantized map and of all 
the Hecke operators "Hecke eigenfunctions" . Our second goal is to 
study the semiclassical limit of the Hecke eigenfunctions. We will 
show that they become equidistributed with respect to Liouville 
measure, that is the expectation values of quantum observables in 
these eigenstates converge to the classical phase-space average of 
the observable. 



I. Introduction 

1.1. Background. One of the key issues of "Quantum Chaos" is the 
nature of the semi-classical limit of eigenstates of classically chaotic 
systems. When the classical system is given by the geodesic flow on a 
compact Riemannian manifold M (or rather, on its co-tangent bundle), 
one can formulate the problem as follows: The quantum Hamiltonian is, 
in suitable units, represented by the positive Laplacian —A on M . To 
measure the distribution of its eigenstates, one starts with a (smooth) 
classical observable, that is a (smooth) function on the unit co-tangent 
bundle S*M, and via some choice of quantization from symbols to 
pseudo-differential operators, forms its quantization Op(/). This a 
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zero-order pseudo-differential operator with principal symbol /. The 
expectation value of 0p(/) in the eigenstate ip is (Op(f)ip, 

Let ipj be a sequence of normalized eigenf unctions: Aipj + = 0, 
Im l^jl 2 = 1- The problem then is to understand the possible limits as 
Xj — ► oo of the distributions 

(1.1) /er^H^,^). 

In the case that the geodesic flow is "chaotic" , it is assumed that the 
eigenfunctions are "random" , for instance in the sense that the expec- 
tation values converge as — > oo to the average of / with respect to 
Liouville measure on S*M. The validity of this for almost all eigen- 
modes if the classical flow is ergodic (so a very weak notion of chaos!) is 
asserted by "Schnirelman's theorem" |2l] [], a fact sometimes referred 
to as quantum ergodicity. The case where there are no exceptional 
subsequences is called Quantum Unique Ergodicity (QUE). Its validity 
seems to be a very difficult problem, to-date unsolved in any case where 
the dynamics are truly chaotic (see however Marklof and Rudnick fl~6|j 
where QUE is proved for an ergodic, though non-mixing, model case). 

1.2. Cat maps. In order to shed some light on the validity of QUE, 
we look at a "toy model" of the situation - the quantization of linear 
hyperbolic automorphisms of the 2-dimensional torus T 2 . Here the 
phase space T 2 is compact and instead of a Hamiltonian flow we con- 
sider a discrete time dynamics, generated by the iterations of a single 
map A G SL(2, Z). If A is hyperbolic, that is \ trA\ > 2, then this 
map is a paradigm of chaotic dynamics. Such maps are sometimes 
called "cat maps" in the physics literature. A quantization of these 
"cat maps" was proposed by Hannay and Berry |J and elaborated on 
in [|T3|, ^, 0, [12], ^Sj . We review this in some detail in Sections 0, [| In 
particular the admissible values of Planck's constant are inverse inte- 
gers h = 1/N, and the Hilbert space of states Hn — L 2 (Z/NZ) of the 
quantum system is finite dimensional, of dimension N = h~ l . To every 
classical observable / G C°°(T 2 ) one associates an operator Op N (f) on 
T~Cn, the corresponding quantum observable. The quantization of the 
cat map is a unitary operator Un(A) on Hn, the quantum propaga- 
tor, unique up to a phase factor, characterized by an exactQ version of 
Egorov's theorem 

(1.2) U N (A)~ 1 Op N (f)U N (A) = Op N (f o A), V/ G C°°(T 2 ) 



1 see Zelditch Q and Colin de Verdiere § for proofs. 

2 This exact version of Egorov's theorem is very special and is a consequence of 
the map being linear. 
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The eigenvectors (f> of the quantum propagator Un(A) are the ana- 
logues of the eigenmodes of the Laplacian, and to study their concen- 
tration properties one forms the distributions 

/ - <0 Pjv (/)^) 

In particular we want to understand the quantum limits as N — ► oo. 
An analogue of Schnirelman's theorem in this setting was proven in 
H One would like to know if QUE holds, that is if the only 

quantum limit is the uniform measure on T 2 . 

The spectrum of the quantum propagator Un(A) has degeneracies, 
which renders the study of possible quantum limits difficult. The de- 
generacies are systematic and are inversely related to the order of A 
mod 2N. Degli Esposti, Graffi and Isola showed that if instead of 
looking at all integer values of N, one restricts to the sparse^ subse- 
quence consisting of primes for which the degeneracies are bounded, and 
moreover split in the quadratic extension of the rationals containing the 
eigenvalues of A, then the only limit is indeed the uniform measure. 

Our first goal in this paper is to show that the degeneracies are 
coupled to the existence of quantum symmetries. There is a commu- 
tative group of unitary operators on TL N which commute with Un(A) 
and therefore act on each eigenspace of Un(A). We will call these 
"Hecke operators" , in analogy with the setting of the modular surface^ 
p0| , 13] , prCjl . We may thus consider eigenfunctions of the desymmetrized 



quantum map, that is eigenstates of both Un(A) and of all the Hecke 
operators. We call these Hecke eigenfunctions. Our second goal is to 
show that these become equidistributed with respect to Liouville mea- 
sure, that is the expectation values of quantum observables in Hecke 
eigenstates converge to the classical phase-space average of the observ- 
able. 

1.3. Results. We turn to a detailed description of our results. We 
first carry out a systematic study of the quantum propagator. We 
define Un(A) so that it only depends on the remainder of A mod 2N 
and satisfies ( |1.2|) . One gets a projective representation A \— > U^{A) 
of the subgroup of "quantizable" elements in the finite modular group 
SL(2,Z/2NZ). We explain (Section |j) that it can be made into an 
ordinary representation if we further restrict to the subgroup T(4, 2N) 



3 It is an open problem to show that there are infinitely many primes where 
the degeneracy is bounded. This is known assuming the Generalized Riemann 
Hypothesis, which in fact guarantees that a positive proportion of the primes satisfy 
the assumption. 

4 A notable difference between our setting and the modular surface is that there 
one expects few, if any, degeneracies. 
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given by g = I mod 4 for N even, g = I mod 2 for iV odd. Thus 
for A,B G r(4,2iV) we have U N {AB) = U N {A)U N {B). Consequently, 
if AB = BA mod 2N then their propagators commute. This is the 
basic principle that we use to form the Hecke operators. 

Fix a hyperbolic matrix A, which we will further assume lies in the 
congruence subgroup 

T(4) = {ge SL(2, Z):g = I mod 4} 

so that its reduction modulo 2N lies in T(4, 2N) for all N. To find 
matrices commuting with A modulo 2N, we use the connection with 
the theory of real quadratic fields (Section ||): If a is an eigenvalue 
of A, form D = Z[a] which is an order in the real quadratic field 
K = Q(a). There is an D-ideal I so that the action of a on / by 
multiplication has A as its matrix in a suitable basis. Thus the action 
of O on / by multiplication gives us an embedding i : O "—>■ Mat 2 (Z), 
and induces a map i : D/2ND — > Mat2(Z/2iVZ). Under this map, 
the images of elements (3 G 0/2NO whose Galois norm is 1 mod 2N 
lie in 5X(2, Z/2iVZ) and commute with A modulo 2N. If we further 
require that (3=1 mod 4D then we get a group of commuting matrices 
i{j3) G r(4, 2N), whose quantum propagators Un(l((3)) commute with 
Un{A) and with each other. These are our Hecke operators. 

Since the Hecke operators commute with Un(A), they act on its 
eigenspaces, and since they commute with each other there is a basis 
of Hn consisting of joint eigenf unctions of Un(A) and the Hecke oper- 
ators, whose elements we call Hecke eigenf unctions. Our main theorem 
is 

Theorem 1. Let A G T(4) be a hyperbolic matrix, and f G C°°(T 2 ) 
a smooth observable. Then for all normalized Hecke eigenf unctions 
4> G Hn ofU]y(A), the expectation values (Op N (f)(f), 4>) converge to the 
phase-space average of f as N — ► oo. Moreover, for all e > we have 

(O Pjv (/)0, 0) = / f(x)dx + Of^N- 1 ' 4 * 6 ), as N - oo 
Jt 2 

Remark 1.1. It is easy to extend Theorem [1] to give similar results for 
matrix elements of Op JV (/). When iV is such that the degeneracies 
in the spectrum of Un(A) are sufficiently small, this implies as in |7j 
that the expectation values of Op N (f) in all eigenstates converge to 
f T 2 f{x)dx. 

Remark 1.2. The exponent of 1/4 in our theorem is certainly not opti- 
mal, and more likely the correct exponent is 1/2. That is the exponent 
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given in where the problem is reduced to one- variable exponen- 
tial sums, which can be estimated using Weil's theorem - the Riemann 
Hypothesis for a curve over a finite field. 

What we in fact show (Theorem is that if fa, i = 1, . . . , N is an 
orthonormal basis of Ti^ consisting of Hecke eigenfunctions then 



N 



(Op N (f)fa,fa) - / f(x)dx 



i 



< iV 



-l+e 



from which we deduce Theorem [I] by taking an orthonormal basis with 
0i = and omitting all but one term on the LHS. If all terms on the 
LHS are of roughly the same size then we would expect this to give the 
exponent 1/2. 

The proof of Theorem [TJ is reduced to a counting problem in Section 
This in turn comes down to counting solutions of the congruence 

Pi - P2 + fa - Pi = mod ND 

in norm-one elements Pi G O/NO. The number of such norm-one 
elements is 0(N 1+e ) (Lemma ||), and since this equation has 3 degrees 
of freedom, the trivial bound of the number of solutions is 0(N 3+e ), 
Ve > 0. To get any result in Theorem [1] we need to show that the 
number of solutions is 0(N 3 ~ S ) for some 5 > 0, that is any saving over 
the trivial bound would do. This is accomplished in Section |7| where 
we show that the number of solutions is 0(N 2+€ ), the optimal bound. 
Acknowledgements: We thank J. Bernstein, D. Kazhdan, J. Keat- 
ing, J. Marklof, F. Mezzadri, P. Sarnak and S. Zelditch for helpful 
discussions concerning various points in the paper. 
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2. Background on quantization of maps 

In this paper we consider the quantization of linear (orientation pre- 
serving) automorphisms of the torus T 2 = R 2 /Z 2 , that is elements of 
the modular group SL(2, Z), which for the most part will be assumed 
hyperbolic (known as "cat maps" in some of the literature). For this 
we first review a procedure (one of several) for quantization of maps. 

The first to quantize the cat map were Hannay and Berry ||. We 
will follow in part an approach by means of representation theory which 
was developed by Knabe and Degli Esposti, Grafli and Isola || 0. 
See also fl2], |], HI for other approaches. 

2.1. The quantization procedure. We start by describing some 
desiderata for a quantization procedure for a symplectic map A of a 
phase space. In the literature it is customary to distinguish two com- 
ponents of the quantization procedure - a kinematic component and a 
dynamical one. 

In the kinematic component one constructs a Hilbert space Tih of 
states of the quantum system^ and an algebra of operators on the 
space - the algebra of quantum observables. Smooth functions / on the 
classical phase space of the system (that is classical observables) are 
mapped to members Op h (f) of this algebra. To make the connection 
with the classical system, it is required that in the limit h — > 0, the 
commutator of the quantization of two observables /, g reproduce the 
quantization of their Poisson bracket {f,g} = J2j §^"§f~ ~ J!f"J^~ : 

(2-1) \ [Op fc (/) f Vh {g)\ - Op,({/, g}) —> 

(we won't specify the sense of convergence). 

The dynamical part of quantization amounts to prescribing a discrete 
time evolution of the algebra of quantum observables, that is a unitary 
map Uh(A) of Hh, which reproduces the classical map A in the limit 
h — > in the sense that: 

(2.2) U h (A)- 1 Ph (f)U h (A) - Op h (/ o A) — — > 

(this is the analogue of Egorov's theorem). 

In our case, the classical phase space is the torus T 2 . The classical 
observables are smooth functions on T 2 . We will find that Planck's 
constant h is restricted to be an inverse integer: h = 1/N, N > 1. 
The state-space Tih will be = L 2 (Z/iVZ). To each observable 
/ G C°°(T 2 ) we will assign, by an analogue of Weyl quantization, an 



operator Op N (f) on Hn so that (|2.1| ) holds, where convergence is in 



stands for Planck's constant. 
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the space of N x N matrices. The dynamics will be given by a linear 
map A G SL(2, Z) so that x = (q ) G T 2 > Ax is a symplectic map 
of the torus. Given an observable / G C°°(T 2 ), the classical evolution 
defined by A is / i— > / o A, where / o A(x) = f(Ax). It turns out that 
for a certain subset of matrices A, there is a unitary map Un(A) on 
L 2 (Z/NZ) so that an exact form of holds: 

U^A)- 1 Op N (f)U N (A) = Op N (f o A), V/ G C7°°(T 2 ) 

This will be our discrete time evolution. 

Below we describe these procedures in detail. 

2.2. Kinematics: The space of states. As the Hilbert space of 
states, we take distributions ip(q) on the line R which are periodic in 
both the position and the momentum representation. As is well known, 
this restricts Planck's constant to take only inverse integer values. We 
review the argument: Recall that the momentum representation of a 
wave-function ip is 

1 r°° 
yh J-oo 

We then require 

ij( q + l) = W(P + 1) = Wp) 

(one may just require that this hold up to a phase). From periodicity 
in the position representation, we get 

^(Q) = ^2c n e(nq) 

neZ 

where 

e(z) := e 2mz 

In the momentum representation, that is applying Th-, we get 

nGZ 

Now in order that Th^{p+ 1) = Fhipijp) we clearly need \ G Z, that is 
for some integer > 1 that 




In that case we also need 

C n+N = c n 

Thus one finds that h = 1/N and the space of states is finite dimen- 
sional, of dimension = 1/h, and consists of periodic point-masses at 
the coordinates q = Q/N, Q G Z. We may then identify TCn with the 
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iV-dimensional vector space L 2 (Z/iVZ), with the inner product ( • , • ) 
defined by 



N 

QmodN 

2.3. Quantizing observables. Next we construct quantum observ- 
ables: For a free particle on the line, we would take as the basic ob- 
servables the position and momentum operators 

#(?) := #(?), P^(<l) ■= 7^(9) 

(h = h/2-K). For our periodic phase space we take the basic observ- 
ables to be e(q) = e 2m ^ and e(p), which correspond to the phase space 
translations 

e(q)ip(q) = e(q)ip(q), e(p)ip(q) = ip(q + h) 
Corresponding to the commutation relation 

we find that 

e(q)e(p) = e~ 2mh e(p)e(q) 

Writing 

h ■= e(p), t 2 := e(q) 
(so that t 2 t\ = e~ 2mh tit 2 ) we put for n = (n 1; n 2 ) G Z 2 

iTrn-i no 

(2.3) T N (n) := e^^t^t™ 1 

Their action on a wave-function ip e L 2 (Z/iVZ) is 



i7rni no 



(2.4) T^nMQ) = eTe(^ Mg + ni) 

These are clearly of period 2N in n: 

T N (n + 2iVm) = T N {n), n,meZ 2 
The adjoint of T^{n) is given by 

(2.5) T N (n)* = T N (-n) 
They also satisfy 

(2.6) T N (m)T N (n) = e^^T N (m + n) 
where 

u;(m, n) = m\n 2 — m 2 ni 
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Now we can finally construct quantum observables: For any smooth 
classical observable / G C°°(T 2 ) with Fourier expansion 

f(x)=J2fne(n-x), x = (?)GT 2 , 
nez 2 

we define its quantization Op N (f) as 

Op N (f) := Yl fM n ) 

The verification of ( |2.1| ) is an easy calculation using ( |2.6| ). 

2.4. The Heisenberg group. We now digress to connect this con- 
struction to the representation theory of a certain Heisenberg group 

H 2 N- 

For vectors x = (x 1 ,x 2 ), y = (2/1, 2/2) define uj(x,y) := x x y 2 - x 2 y x . 
This is a non-degenerate symplectic form. The Heisenberg group H 2 n 
is defined to be the set (Z/2NZ) 2 x Z/2NZ with multiplication 

(x, z) ■ (x', z') := (x + x', z + z' + lu(x, x')) 

This is at odds with the standard convention where one multiplies u 
by 1/2, but is essential for us because 2 is not invertible in Z/2NZ. 

It is useful to record various facts about the multiplication in H 2N : 
The inverse of (x, z) is 

(2.7) [x,z)- 1 = (-x,-z) 
The commutator of two elements is given by 

(2.8) (x, z)(x', z')(x, z)~ l (x\ z'y 1 = (0, 2u(x, x 1 )) 

From this commutator identity and the fact that uj is non-degenerate 
we immediately find 

Lemma 2. The center of H 2N is (NZ/2NZ) 2 x Z/2NZ, that is 
Cent(H 2N ) = {(Ne, Nrj, z) : e, 77 = 0, 1, zE Z/2NZ} 
We define a representation of H 2 n on L 2 (Z/NZ) by setting 

ir(n,z) = e(—)T N {n) 

From the relation ( |2.6| ) it follows that 7r(h)7r(h') = Tr(hh'), i.e. we do 
indeed get a representation. 

The center of H 2 ^ then acts via the character x given by 

z + x y 

X{x ,y ,z) = e{ 2N ) 

(that is vr(x , y , z) = x(x , y , z)I). 
The basic facts about 7r and the representation theory of H 2N are 



10 PAR KURLBERG AND ZEEV RUDNICK 

Proposition 3. i) All irreducible representations of H 2 n have dimen- 
sion at most N . 

ii) The representation tt is irreducible, and is the unique irreducible 
N -dimensional representation with central character \. 

We omit the details of the proof; the main point (which is easy to 
verify from the definitions) is 



Lemma 4. the trace ofT^{n) is given by 
I tiT N (n)\ = 



N if n = (0,0) mod iV, 
otherwise. 

Proof. Let fa = \fNb~i where 5{ is the Dirac delta function supported 
at i, so that {<pi}f =l is an orthonormal basis of L 2 (Z/NZ). Then 

N 

tr T N (n) = ^2(T N (n)fa,fa) 



i=l 



and by equation ( |2.4j ) 



TNinjfaiQ) = e( — )fa(Q + n h 

f n 1 n 2 + 2n 2 Q. 

-jfa-miQ) 



2N 

—n\n 2 + 2n 2 i 



2N 

Therefore trTjy(n) = unless fix = mod N, in which case 

N N 
i=l i=l 

The result now follows since ^2^ =1 e(^) equals iV if n 2 = mod N, 
and is zero otherwise. □ 

2.5. Description of 7r as an induced representation. Let Y be 

the subgroup of elements 

Y = {(x , y,z) :y,ze Z/2NZ, x G NZ/2NZ} 

It is easily seen to be a normal, maximal abelian subgroup, of index 
N, containing the center. Set for (x , y, z) G Y 

z + x y 

T{x ,y,z) := e( 2N ) 

This is a character of Y (we need to use 2xq = mod 2N in verifying 
this), restricting to the character x{ x o,Vo, z ) — e( z+ 2 ^ yo ) of the center. 
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We consider the induced representation Indy 2JV r of the Heisenberg 
group. The basic model for it is the space of functions $ : H 2 n C 
satisfying $(a/i) = r(a)$(/i) for a G Y , h G H 2 n- The action of the 
group is by right multiplication h$(h f ) := $(h'h). By restricting to the 
subgroup X — {(x, 0, 0)} we can realize this induced representation as 
functions on Z/2iVZ which are iV-periodic (since the element (N, 0,0) 
lies in Xf]Y). This space of functions we can identify with L 2 (Z/iVZ). 

Let us compute the action of a group element h = (x, y, z) G H 2 n in 
this model. For this we need to write (x 1 , 0, 0) • h as a ■ (x", 0, 0), a G Y. 
The relevant identity is 

(x',0, 0) (x, y, z) = (0, y, z + xy + 2x'y) (x' + x, 0, 0) 

Thus the element h = (x, y, z) acts as 

^ ) = e ( 2JV X > 

In particular (x, 0, 0) acts as translation by x and (0, y, 0) as a multipli- 
cation operator 4>{x') \— > e(^)0(x'). The center acts by the character 
(xo,yo, -2) ^ e ( 2 "27v W )- These show that 7r coincides with the induced 
representation Indy 2iV r. 
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3. Dynamics: quantized cat maps 

We now show how to assign to (certain) linear automorphisms A of 
the torus T 2 a unitary operator U N (A) on L 2 (Z/NZ), which satisfies: 
For all observables / G C7°°(T 2 ) 

U^A)- 1 Op N (f)U N (A) = Op N (f o A), 

The finite modular group SL(2, Z/2iVZ) acts by automorphisms on 
the Heisenberg groups H 2 n via (x, z) A := (xA, z), A G SL(2, Z/2NZ). 
That this is indeed an automorphism (that is (hih 2 ) A — h A h 2 ) follows 
from A preserving the symplectic form uj. Moreover we have (h A ) B = 
h AB . Composing the representation it of H 2 n with A gives a new 
representation n A (h) := ir(h A ), which is clearly still an irreducible 
iV-dimensional representation. Its central character x A can be easily 
computed as follows: if x , y G NZ/2NZ and (xi, y-y) = (x , yo)A then 
X A is given by 

X A (x , y , z) = x((x , yo)A, z) = e( - ^ Vl ) 

This will be the same character as x iff x iUi = x oVo mod 2N for 
all x ,y G NZ/2NZ. Writing A — ( a c b d ) and x — Ne, y — Nrj, 
e,i] G Z/2Z, this is equivalent to requiring 

N{abe 2 + cdrf) = mod 2, Ve, r] G Z/2Z 

or 

Nab = Ned = mod 2 

This is only a restriction if N is odd, and is satisfied by the elements 
of the theta group 

r e (2iV) = {(l ^ e SL(2, Z/2NZ) :ab = cd = mod 2} 

Therefore if A G r (2A^), we get a unitarily equivalent representa- 
tion n A of H 2 n ■ Thus there is a unitary map Un(A), the quantum 
propagator associated to A, so that 

tt(/i a ) = U N (Ay\(h)U N (A), V/i G fT 2JV 

In particular we find 

(3.1) ^(^-^(rOE^A) = 7V(nA) 
and consequently for all observables / G C°°(T 2 ), 

(3.2) Op N (f o A) = U N (A)~ 1 Op N (f)U N (A) 

We now for any "quantizable" element A G SL(2, Z) (that is A = 
(cd) w ith ab = cd = mod 2), we define the quantum propagator 
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(or "quantized cat map") to be U N (A) where A e SL(2, Z/2NZ) is 
the reduction of A modulo 2N. Thus by its construction, Un{A) only 
depends on the reduction A mod 2N. (This is a difference from the 
construction in Hannay and Berry ||). 
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4. MULTIPLICATIVITY 

The quantum propagators Un(A) are uniquely defined up to a phase- 
factor, because of the irreducibility of tt (Schur's lemma). Thus they 
define a projective representation of T 9 (2N), that is 

U N (AB) = e t4>N{A > B) U N (A)U N (B) A, Be T e (2N) 

Define the subgroup 

r(4, 27V) = L e SL(2, Z/2NZ) : ( 9 = ' *»* *■ " ™) 

[g — I mod 2, N odd I 

The goal of this section is to show that there is a choice of phases 
for the propagators Un(A) so that on the subgroup r(4, 2N) the map 
A i — > Un{A) is a homomorphism: 

Theorem 5. There is a choice of quantum propagators so that 

U N (AB) = U N (A)U N (B) } A, Be T(4, 2N) 

As a consequence we find 

Corollary 6. If A,B 6 r(4, 2iV) commute mod 2N then their propa- 
gators also commute: Un(A)Un(B) = Un(B)Un(A). 

Theorem [5] in various guises is essentially known, and arose out of 
the study of theta-functions and the Weil representation. One form is 
due to Kubota [Q (see also ||). There are also treatments purely at 
the finite level 0, Since Corollary |] is absolutely crucial to our 
work, and we did not find a good reference for the exact form that 
we need, we will sketch a proof (or more precisely, a verification) of 
Theorem [5[ We wish to note that Theorem |5| is a-priori more subtle 
than Corollary since once we know that there is some choice of phases 
for which Corollary |6| holds, than it holds for all choices; this is not the 
case with Theorem [|.(] 

4.1. Reduction to prime powers. Factor 2N = Y[ p P kp = 2 fc Yl P >2P kp = 
2 k M, M odd. The Chinese remainder theorem gives an isomorphism 

Z/2NZ ~Y[Z/p kp Z 

v 

given by 

x i — ► (x mod p kp ) p 



6 We thank Jon Keating for emphasizing this point to us. 
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with inverse 



mod p kp ) p i— > —jrr p Xp mod 2N 



where r p is the inverse of 2N/p kp modulo p kp . Correspondingly we have 
a bijection 

L 2 {Z/2NZ) ~ (g)L 2 (Z/p kp Z) 

p 

We define the phase space translations on L 2 {Z/p kp Z) as in ( |2.4| ), 
by 

rW W (Q) = e( r ' ( "'"% +2 "^ ) W + ni ) 

Is is then a simple matter to see that T^{n) = ® p T^(n), that is if 
ip = ® p 4>p G (g) p L 2 (Z/p kp Z) is decomposable then 

V 

This allows us to express the quantum propagators Un(A) as tensor 
products. Indeed, if we already have propagators (A) which satisfy 

(4.1) U {p \A)- l T^\n)U {p \A) = T {p \nA) 
We then set 

(4.2) U N (A) := 0U (P \A) 
which still satisfies 

U N {A)- l T N (n)U N {A) = T N (nA) 

for all n G Z 2 and therefore Un(A) coincides up to a phase with any 
other map satisfying this. 

We use this procedure to define Un(A) (that is, choose a phase) 
so that Un is an honest representation of a subgroup r(4, 2N) of 
SL(2, Z/2NZ), not merely a projective representation. From the fac- 
torization property (fO|), it follows that it enough to show that U {p) is 
a representation of SL(2, Z/p kp Z) when p > 2 is odd, and of T(4, 2 k ) if 
N = 2 k ~ 1 M is even. 

4.2. Gauss sums. We need some preliminary information on Gauss 
sums. We define normalized Gauss sums 

(4.3) S P (a,p*) = J= Yl < rVar2 



P imodp* ' 
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For p odd these are 4-th roots of unity. To describe them, define for 

t e (Z/ P k zy 

_ S r (t,p k ) 

Ar * h[t) - S r (l,p k ) 

Note that if t = t\ e (Z/p k Z)* is a square then A rp k(t) = 1 since from 
(|4.3| ) we find after the change of variables x± = t±x that S r (t,p k ) = 

S r (l, P k ). 

For odd p, A rp k is given in terms of the Legendre symbol as 

/t\ k 

tV^, pk \ ij ~ 



and is a character of {Z/p k Z)*\ 
When p = 2, one has 

KAt) = ( z f) *- r{p - 1)/8 

where t is the smallest positive residue of t mod 4. In that case it is 
not quite a character of the whole multiplicative group of Z/2 k Z, but 
instead satisfies 

(4.4) A r . 2k (tf) = (t, t') 2 A r . 2 k (t)A r<2k (£') 

where (t, t') 2 is the Hilbert symbol. In particular, if t, t' — 1 mod 4 then 
the Hilbert symbol is trivial and so we get a character of the subgroup 
{t = 1 mod 4} C (Z/2 fc Z)* [] given simply by 



A r2 k(t) 



1, t = 1 mod 8 

(-l) k , t = 5 mod 8 



For p odd we will also need to know the normalized Gauss sum ( |4.3| 
when t — — 1 in which case one has 



S r (-l,p k ) 

where 

e(p) = 



1, k even 

e(p)(V), k odd 

1, p = 1 mod 4 
i, p = 3 mod 4 



7 this is relevant for k > 2. 
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4.3. p odd. We describe how to define on SL(2, Z/p k Z) so that 
it gives a representation - see Nobs |l8j for details. This group is 
generated by the matrices 

1 % f ( A L 1 



(4.5) 



and so it suffices to specify on such matrices, provided we preserve 
all relations between them. This is done by the formulas 

(4.6) l/wf 1 fW) = e(^MaO 

(4.7) U^f* t _ i y )l>(x) = A T A t M tx ) 



2rxy, 

el 



(4.8) I/Wf x 1 ^(x) = S r (-l,p k )^= J2 MvU pk 

vP ymodp k 

It is easy to check that these satisfy ( |4.1[ ). To see a verification that 
this prescription does indeed give a consistent definition (that is that all 
relations between the generators (|4.5f) are satisfied), see e.g. jl8f . Once 
we have this then automatically we get U ip) (AB) = U ip) (A)U {p) (B) . 

Remark 4.1. It is in fact the case that any projective representation 
of SL(2, Z/p k Z), p odd, can be modified to give a representation (and 
more generally, SL(2, Z/mZ) if m ^ mod 4)- this is due to Schur f22|j 
when k — 1. See |L7] and [0 for the general case. 

4.4. p = 2. Here we restrict to the subgroup T(4, 2 fc ), k > 2. The 
literature in this case is harder to come by, so we include complete 
proofs. We start by describing generators and relations for this group. 
More generally, let p be any prime, and k > 2. Let 

Y{p\p k ) := {g E SL(2, Z/p k Z) : g = I mod p 2 }. 

Lemma 7. T(p 2 ,p k ) has a presentation with generators u + (x), U-(y), 
s(t), where x,y,t G Z/p k Z, x, y = mod p 2 , t = 1 modp 2 , and rela- 
tions 

(4.9) u + (x)u + (x') = u + (x + x 1 ) 

(4.10) u_(y)u_(y') = u^y + y') 

(4.11) s(t)s{t') = s(tt') 

(4.12) sitfu+Wsii)- 1 = u + (t 2 x) 

(4.13) s{t)u^y)s{ty l = u-(r 2 y) 

(4.14) s(d)u + (a)u-(b) = M_(^ 1 6)w + (rfa), d := (1 + ab)~ l 
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Proof. Let G be the abstract group with the above presentation. We 
get a map ^ from G into T(p 2 ,p k ) by taking 

* : u+(x) fj, u-{y) ~Q A, s(t) ^ P rl 

One verifies that the relations hold in SL(2,Z/p k Z) so that $ is a 
homomorphism. Next, note that we have a "Bruhat decomposition" 
for T(p 2 ,p h ): Every element can be uniquely written in the form 

't \ (\ x\ (\ 
which follows from the formula 



a b\ id \ 1 1 bd 



c d) 1 \ d) \ 1 

(note that since d — 1 mod p 2 it is in particular invertible). This 
implies that the map \I/ is surjective. To see that \l/ is an isomorphism, 
it suffices to show that every element of the abstract group G can also be 
written in the form g = s(t)u+(x)u-(y) , since then by the uniqueness 
of the decomposition in T(p 2 ,p k ), \l/ is also one-to-one. 

With the aid of the first five relations, every word W G G can be 
written product: 

W = s(ti)w + (a;i)«_()/i) s(t n )u + (x n )u-(y n ) 

for some n > 1. We prove by induction on n that we can write W = 
s(t)u + (x)u-(y) for x,y = mod p 2 , t = 1 modp 2 . When n = 1 this 
holds trivially, and for n > 1 we use the relation (|4.13| ), ( |4.14|) to write 

u-(y n -x)s(t n )u + (x n ) = s(t n )M_(t 2 ?/ n _i)M + (x n ) = s(t n )s(t')u + (x')u-(y') 

and so 

W = s(ti)u + (x 1 )n_(?/ 1 ) . . . s(t n _i)n + (s n _ 1 )s(t n )s(t / )«+(^ / ) w -(z/ / ) M -(z/n) 

= s(ti)u+(xi)u_(3/i) . . . s(t' n _ 1 )u+(x"_ Ou-G/^i) 

after a further application of the first five relations. The result now 
follows by induction. □ 

We now specify the propagators U <y2 \A) for the generators: For 

and ^ they are given by the same formulas Q4.6p, ([4.71). 

For the matrices 
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we conjugate ( |4.6| ) by an analogue of the Fourier transform ( |4.8|) and 

define 

(4.15) 

ymod2 fe 2mod2* 

To show that this defines a representation, one has to check that all 
the relations of Lemma [7] are satisfied. The first five are fairly straight- 
forward, bearing in mind that A is a character of the multiplicative 
group of residues t = 1 mod 4 (see ( }4.4| )). The last relation ( f4.14| ) 
requires verifying an identity of Gauss sums: Unwinding the action of 
the right and left hand sides in fl4.14| ) we must show that 



A (d) ^(y)e^(2yz-bz 2 - 2dxz + ad 2 

z mod 2 k y mod 2 k 

= E E ^(y)e(^{2yz-d~ 1 bz 2 ~2xz + ady 2 )^ 

z mod 2 k y mod 2 k 

Now, d = 1 mod 16 implies that A(d) = 1 since then d is a square 
modulo 2 fc , and if the identity is to hold for all ip and all values of x 
we obtain that for all x, y 



(4. i6) e (^k (~ bz2 + 2z (y - dx "> + ad2x2 ) ' 

z mod 2 k 

= £ e (- rf ~ lfe2 + 2 «(f - x ) + ac ^ 2 ) ) • 

a mod 2 k 

We will verify this in Appendix |A|. 
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5. Hecke operators 

We now introduce a commutative group of unitary operators on 
L 2 (Z/iVZ) which commute with Un(A). For this, we have to bring 
in the theory of quadratic fields (see [HJ for a survey in connection to 
cat maps). 

5.1. Integral matrices and quadratic fields. Let A G SX 2 (Z) be a 
hyperbolic matrix: \trA\ > 2. The eigenvalues a, a" 1 of A generate a 
field extension K = Q(a), which is a real quadratic field since tr(A) 2 > 
4. We denote by Ok the ring of integers of K. The eigenvalues a, a" 1 of 
A will be units in Ok- Adjoining a to Z gives an order O = Z[ot] C iD^ 
in X. We claim that there is an O-ideal / C D so that the action of 
a by multiplication on / is equivalent to the action of A on Z 2 , in the 
sense that there is a basis of / with respect to which the matrix of a 
is precisely A. 

The construction is as follows [p3|| : Since a is an eigenvalue of A, 
there is a vector v = (t>i,t>2) such that vA = av and v G D 2 . Let 
/ := Z[v i, t>2] C £). Then J is in an D-ideal, and the matrix of a acting 
on I by multiplication in the basis t>i,t>2 is precisely A. 

Remark 5.1. It is easy to check that the above construction sets up 
a bijection between GL 2 (Z)-conjugacy classes of elements in SL 2 (Z) 
with eigenvalues a, a -1 and and ideal classes in the order D. (Recall 
that two ideals I\, J 2 are said to be in the same ideal class if there exist 
nonzero a, b G O so that al\ = 6/2-) 

In the same way, the action of O by multiplication on / gives us an 
embedding 

l: D Mat 2 (Z) 

so that 7 = x + ya G corresponds to xl + yA. Moreover, the 
determinant of xl + yA equals J\f(j) = 77, where Af : K — > Q is the 
Galois norm. In particular, if 7 G D has norm one then 7 corresponds 
to an element in SX 2 (Z), and if in addition 7=1 mod 4D then 7 
corresponds to an element in T(4). 

5.2. Hecke operators. Given an integer M > 1, the embedding t : 
D Mat 2 (Z) induces a map t M : C/MD Mat 2 (Z/MZ) and the 
norm M : K — > Q gives a well-defined map 

M-.O/MD -> Z/MZ. 

We let C^(M) be the group of norm one elements in O/MO. 

C A {M) = ker [JV : (O/MO)* -> (Z/MZ)*] 
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Similarly, replacing the order O by the maximal order Ok we set 
C K (M) = ker [AT : (D K /MD K )* (Z/MZ)*] 

to be the norm one elements in Ok/MOk- 

If M = 2N is even, we set C e A (M) to be the elements of C A (2N) that 
are congruent to one modulo 4D (respectively 2D) if N is even (resp. 
odd). For M odd we set C e A {M) = C A (M). 

By construction, the image of C e A (2N) in Mat 2 (Z/2iVZ) lies in T(4, 2N). 
Since a commutes with all elements in C A (2N) we see that A commutes, 
modulo 2N, with the elements in l(C a (2N)) . Thus by Corollary |6| the 
quantizations Un(l(/3)) of (3 G C A {2N) commute with Un(A) and with 
each other. We will call these "Hecke operators" . 

We will need to know the number of Hecke operators: 

Lemma 8. The number of elements of C A (2N) satisfies 
N 1 -' < \C A {2N)\ < N 1+t , Ve > 0. 

Proof. Since the reduction map D — > D/4D has image of size 4 2 , 
C A (2N) has bounded index in C A (2N). The inclusion O C Ok in- 
duces a map O/MD — > Dk/MDk which has kernel and co- kernel 
of size at most [Ok '■ £) ], independent of M. Therefore the induced 
map C A (M) — > Ck{M) on norm-one elements also has bounded kernel 
and co-kernel. Thus it suffices to prove the lemma in the case of the 
maximal order Ok- By the Chinese remainder theorem, it suffices to 
prove it in the case of prime powers, which is given in Appendix p] by 



Lemma 19. □ 



5.3. Hecke eigen-functions. The Hecke operators Un(i((3)), (5 G 

C e A (2N), commute with each other and with Un(A). Therefore the 
eigen-spaces of the unitary map Un(A) break up into joint eigen-spaces 
of the Hecke operators. Such a joint eigen-f unction we call a Hecke 
eig en- function. In other words, there exist an orthonormal basis {(f>i} 
of L 2 (Z/iVZ) and characters Aj of C A (2N) such that 0, are eigenfunc- 
tions of Un(A) and 

UnW))^ = \(0)4>i, V/3 G C e A (2N). 

We call such a basis of L 2 (Z/NZ) a Hecke basis. 
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6. Ergodicity OF Hecke EIGENFUNCTIONS 

In this section and the following we show that if G L 2 (Z/NZ) is a 
normalized Hecke eigenfunction then the expectation values (Op Ar (/)0, 4 
converge to the classical phase-space average J T2 / for all smooth ob- 
servables (Theorem [I] of the Introduction). In fact, we show something 
stronger: 

Theorem 9. Let fa G L 2 (Z/NZ) , i = 1,...,N be any orthonormal 
basis of Hecke eigenf unctions ofU^{A). Then 



N 

E 

i=l 



(Op N (f)fa,fa) - I f(x)dx 

T 2 



-1+6 



6.1. Proof of Theorem To prove this theorem, it suffices (see 
below) to prove it for the basic observables f(x) = e(nx), ^ n G Z 2 , 
that is to show 

Theorem 10. Let 0/n6 Z 2 ; and let fa G L 2 (Z/NZ), i = 1, . . . , N 
be any orthonormal basis of Hecke eigenf unctions ofUN{A). Then 

N 

J2\(T N (n)fa,fa)\ 4 <t: e \n\ ie N~ 1+ % N -> 00 
i=i 



The proof of Theorem ^| from Theorem 10 is easy using the rapid de- 
cay of the Fourier coefficients of /. Indeed, write f(x) = J2nez 2 f( n ) e { nx )i 
so that Op N (f) = Ytnez? K n ) T N(n). Therefore 



N 

E 



4 N 



E 

i=l 



f( n )( T N(n)fa,fa) 



o^nez 2 



N 



<Y1 IIi^ nfe )^ Tjv ( nfe )^'^) 

i=l n± ,...,n4^0 k=l 



For notational convenience we write 



ti(n) := \(T N (n)fa,fa)\ 
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Now interchange the order of summation, and apply Cauchy-Schwartz 
twice: For fixed ni, 712, n,3, n.4 



N 



< 



i=i 

N 



V2 / N v 1/2 4 , N 



1/4 



< 



n e^w 



i=l 



i=l 



k=l \i=l 



Now use Theorem [10]: For nj; ^ 0, 

N \ V4 



5> 



,i=i 



and so we get 



N 



i=l 



fc=l 



Now sum over all possible nk 7^ to find 



N 

£ 

i=l 



(Ov N (f)(f>i,(f>i) - / /(aOcfo 



<iV- 1+£ (^/(n)H 



which proves Theorem 



□ 



6.2. Reduction to a counting problem. We first reduce Theo- 
rem 



10 to a counting problem. 



Proposition 11. Fix 7^ n = l(u) e Z 2 , v £ I . Then for any Heche 
basis of eigenf unctions (f>i, 



N 



Y,\{TN{n)4>iAi)\ A < 



i=i 



N 



e C e A {2N) : v{Px - ft + (h - (3 A ) = mod NI} 



\C e A (2N)\* 

In order to prove Proposition [11], we define for n = 7^ v G I 

1 



D = D(n) 



|C*(2iV)| 



/3eC A (2N) 



24 PAR KURLBERG AND ZEEV RUDNICK 

If (tij) is the matrix coefficients of Tjy(n) expressed in the eigenvector 
basis {4>k\ so that Uj =< Tjv(n)0j, <pj >, then we see that 



Since the sum of a nontrivial character over all elements in a group 



/36C«(2n) 



vanishes we have 

(6.1) D i:j 



tij if Aj — Aj, 
otherwise. 



Lemma 12. Wii/i D defined as above we have 

£ |%| 4 <tr((M>) 2 ). 

Proof. Let £) = (djj) = (vj) where the f,'s are the column vectors of 
D. Examining the (k, fc)-entry of (D*D) 2 we get 

{{D*D) 2 ) kk = < v u v k >< v k , Vi >= | < fj, f fc > | 2 , 

i i 

and hence 

tr((£>*D) a ) >Y,\< v k,v k >\ 2 >J2 Idyl* 

The result now follows from equation (|6.1| ). □ 
Lemma 13. VFe have 

tr((D*D) 2 ) < lc e^ N){4 1 {A G C*(2iV) : ^(ft-ft+ft-ft) = mod iV/}| 



Proof. Recall that by (pJj) , since n • = i(z//5) for /5 G D, n = t(u), 

U N (i(f3)r 1 T N (n)U N (i((3))=T N (i(uP)) 

Also note that T N (w)* = T N (—w) for all w by ( |2.5| ). Substituting 
the definition of D and expanding we see that (D*D) 2 is given by 
l/\C e A (2N)\ 4 times a sum, ranging over all ft, ft, ft, ft G C^(2iV), of 
terms 

T J v(i(^ 1 ))r J v(- t (^ 2 ))T J v(^/3 3 ))Tiv(-^4)) 

= 7(A, ft, ft, ft)2V(i(i/(ft - ft + ft - ft))) 

where 7(ft,ft,ft,ft) has absolute value one (see ( |2.6| )). Now take the 
trace; by Lemma |, the absolute value of the trace of T^{n) equals N 
if n = (0, 0) mod N, zero otherwise. The result now follows by taking 
absolute values and summing over all ft, ft, ft, ft G C^(2iV). □ 
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It remains to estimate the number of solutions of 
(6.2) + =0 mod NI, ^eC° A {2N). 

We will show 



Proposition 14. The number of solutions to equation ( |6.2| ) is bounded 
byO{\Ar{v)\ 8 N 2+ <). 



we 



6.3. Proof of Theorem |10f Conclusion. By Proposition [XT 
need a suitable upper bound for the number of solutions of equation 
and a lower bound for the number of elements of C A {2N). By 



Proposition 14, the number of solutions is at most \M N 2+e . Note 
that \N{v)\ < \n\ 2 . From Lemma || we obtain that \C e A {2N)\ > N 1 ~ e 
and the result follows. 
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7. Counting solutions 



In this section, we prove Proposition 14. 

7.1. A reduction. Since NI C iVD C NOk, the number of solutions 
to ( p.2|) is bounded by the number of solutions to 

Moreover, at the cost of increasing slightly the number of solutions, we 
may omit the parity condition on /3j and so replace C e A {2N) by Ca(2N). 

The inclusion O C Ok induces a map O/MO — > Ok/MOk which 
has kernel and co-kernel of size at most [Ok '■ O], independent of M. 
Therefore the induced map 

C A {M) = ker[(D/MD)* -> (Z/MZ)*] 

-> C K (M) = ker[(D*/MD x )* - (Z/MZ)*] 

on norm-one elements also has bounded kernel and co-kernel. Thus, up 
to a bounded factor (depending on A but not on iV or u), the number 
of solutions to ( |6.2| ) is bounded by the number of solutions of 

(7.1) v(Pi-P2+fo-P*) = mod NO K , Pi G C K (2N) 

At the cost of increasing the number of solutions, we multiply the 
equation (|7.1|) by the Galois conjugate v to get an equation 

M(v)(j3 1 -/3 2 + fo-Pi) = mod NO K , Pi G C K (2N) 
Setting 



gcd(i\r,jV(>/)) 

this equation is equivalent to 

(7.2) /3i-/3 2 + /3 3 -& = modiVD^, PieC K (2N) 

Next, note that the reduction map Ok/tsOk — > Ok/tOk has ker- 
nel rOK/rsOK — Ok/sOk of size s 2 , and so the induced map on 
norm-one elements C#-(rs) — > C^(r) has kernel of order at most s 2 
(this is crude, but sufficient for our purposes). Thus the reduction 
map Ck{2N) -> C^iV') has kernel of size at most 4gcd(iV, N(v)) 2 < 
A\N(v)\ 2 . Therefore the number of solutions of ( |7.2| ) is bounded by 
(4|A/"(z/)| 2 ) 4 times the number of solutions of the equation 

(7.3) Pi-P2 + P3~Pt = mod N'O k , Pi E C k (N') 



Equation (7j3) is invariant under Galois conjugation and we obtain 
a second equation (note that j3 = P' 1 since N(0) = 1 mod N') 

(7.4) flf 1 - fe 1 + Pz 1 -P^ = mod N'O k . 
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7.2. A transformation. We thus have a system of equations ( ]7.3| ). 
(|7.4| ), which we transform using the following: 

Lemma 15. Ifx,y,z,w are invertible then the system of equations 

x + y = z + w 
x^ 1 + y^ 1 = z^ 1 + w~ l 

is equivalent to the system 

(z- x)(z-y)(x + y) = 
w = x + y — z 

Proof. From the second equation we get 

x + y z + w 
xy zw 

or 

(x + y)zw — (z + w)xy 

The first equation gives that w = x + y — z, inserting it in (x + y)zw = 
(z + w)xy we get 

(x + y)z(x + y — z) = {x + y)xy 

or 

= (x + y)(zx + zy — z 2 — xy) = — (z — x)(z — y)(x + y). 

□ 



Thus by lemma [15] the system of equations (|7.3|) , (|7.4j ) is equivalent 
to the system: 

(7.5) (ft-0i)(/33-#O(/3i+#O=O modiV'D^ 

(7.6) ft = ft - (3 2 + p 3 mod N'Q K . 

with ft G C^(iV'). 

Since ft is determined by ft, ft, ft, we may ignore the second equa- 
tion ( [7.61) (at the cost of increasing the number of solutions, since being 
in Ck(N') is a non-empty condition). Multiplying equation (|7.5| ) by ft" 3 
and letting ft' = ft /ft we obtain 

(7.7) (i_ft)(i_ft)(ft+ft) = modiV'D^ 



Since ft is arbitrary, the number of solutions of ( [7.5|) is bounded by 
|Cjr(iV')| times the number of solutions in ft, ft G C K (N') to ( ff7|) . 
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7.3. Prime powers. By the Chinese remainder theorem the number 
of solutions to (|7.7|) is multiplicative, and we may concentrate on the 
prime power case. Thus we need to count the solutions to the equation 

(7.8) (l-#)(l-#)(# + $) = modp k O K 

with (3[ G K /p k O K , U{(3[) = 1 mod p k . 

We first recall some properties of primes in quadratic extensions: Let 
P\p be a prime in Ok lying above p, and let e denote the ramification 
index, i.e. the largest integer e such that P e \pO K . Since K is quadratic 
e G {1,2}, and e = 1 for all but finitely many primes p. If e = 2 then p 
is said to be ramified. If e = 1 then p is called unramified, and one of 
two things can happen: either p£>K — P is still a prime ideal, in which 
case p is said to be inert, or pOx = PP, in which case p is said to split. 

Now, fix a prime p with ramification index e, be it one or two. The 
norm map Af : Ok — > Z gives a well-defined homomorphism 

(o K /p ek ) x -> (z//) x . 

We let 

(d^/p 6 *) 1 

be the kernel of this map, i.e., the group of norm one elements. For 
I < ek we let 

((l + p'ytl + F*)) 1 

be the norm one elements in the subgroup (1 + P l ) /(l + P efc ), these are 
precisely the norm one elements that reduce to one modulo P l . 

Lemma 16. There is a constant c> 1 so that the number of solutions 
of equation (|7.8j ) is at most ckp k . 

Proof. Equation ( |7.8| ) is invariant under Galois conjugation, therefore 

its solutions in Ox/p k OK correspond bijectively to solutions (3[ G 

Ox/P ek , N{(3'i) = 1 mod p k (this is of course only an issue in the 

k 

split case where Ok/p 'Ok — Ok/P x Ok/P )• Thus we need to 
count solutions of 

(7.9) (l-^)(l-^)(^ + ^) = modP efc 

with (3[ G K /P ek , JV(#) = 1 mod 

We will first assume that p is o<i<i Since (3[ = f3' 2 = 1 mod P implies 
that /3[ + /?2 = 2 ^ mod P we see that at most two of the factors in 
equation ( |7.9|) can be congruent to zero modulo P. Moreover, we may 
assume that the third factor is nonzero by multiplying by a suitable 
(3 and permuting the variables. (Of course we must then compensate 
by multiplying the number of solutions by fj!) ) . Now, if the product is 
zero modulo P ek , then there is some < n < ek such that one factor is 
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zero modulo P n , and the other zero modulo P ek n . Thus the number 
of solutions to equation ( [7.9|) equals 

efc-1 

. _ . _ , .11 I. ... , .. _ , . Ll + 

n=l 



|((! + P n )/(1 + P efc )) x ((l + P efc - n )/(l + P efe )) 

(^/p^) 1 



Using Lemma we obtain 

((1 + P n )/(1 + P efe )) 1 x ((l + P efc - n )/(l + P efc )) 1 
and by Lemma |i~9| 

{0 K /P ek ) 1 <2(p + l)p k -\ 



< 



Hence for p odd, the total number of solutions to (|7 
4(p + + 3(e£; - l)p e -y < 

1 for all but finitely many primes). 



is bounded by 



{ since e 

If p = 2 it is no longer true that only two factors can be zero mod- 
ulo P. However, /3i = /3 2 = 1 mod P e+1 implies that /?i + fa = 2 
mod P e+1 . Since 2D^ = P e , we see that if two factors are zero modulo 
P e+1 , then the third factor can be congruent to zero at most modulo 
P e . We may thus bound the number of solutions by counting the num- 
ber of ways the product of two factors can be equal to zero modulo 
pefc-e T;hj s we can do as we f or odd primes, and we obtain the 
same bound as before, except that we lose an additional factor of at 
most 



((1 + P ek - e )/{l + P ek )Y <2° (e) 
This proves the Lemma. 



Oil) 



□ 



7.4. Proof of Proposition [14] . By multiplying over all primes, we 
see from Lemma [16] that the number of solutions of equation ( |7.7| ) is 
0((N') 1+e ). Therefore we see that the number of solutions of ( |7.5| ) 
is 0((N') 2+t ) since \C K (N')\ < (N') 1+e by Lemma 0. This gives a 
bound for the solutions of (|7.3| ) and multiplying by |A/"(z/)| 8 gives a 
bound for the number of solutions of ( |7.2| ). In turn, by the reasoning 
in section [TTTJthis gives a bound of 0(|A/"(z/)| 8 iV 2+<: ) on the solutions of 
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Appendix A. An identity of Gauss sums 

For section |] we need to prove the identity ( 4.16 ). To prove it we 
will need a lemma about Gauss sums. Given an integer x, we define its 
"dyadic valuation", v(x), by x = 2 v ^x , where x is an odd integer. 
Let 

G(b,c)= e(^{-bz 2 + 2cz)). 

z mod 2 k 

Lemma 17. If v(c) < v(b) < k then 

2 k if v (b) = k-l andv(c) = k-2, 



otherwise. 



G{b, c) 
Proof. We may write 

G(b,c)= Yl e (|r(-^ 2 + ^)) 



mod 2 k 



where (3 is an integer satisfying 2c/3 = b mod 2 . Let n = k — l — v(c); 
it is the smallest integer n such that e (%pr%) = 1 for all x = mod 2 n . 

Assume first that n > 1. Let e = e 2™~ 1 be such that e (t§tc) 7^ 1. 
Making the change of variables z — > z + e we see that 

G(b,c)= e(^(-P(z 2 + 2ez + e 2 ) + z + e)) 



2 mod 2 fc 

'2cr 
"2 



G(&,c)e(^e 



since 2ez + e 2 = mod 2 n . But e (f£e) ^ 1 and therefore G(6, c) = 0. 

If n < 1 then, as n = k — 1 — v(c) and v(c) < t> (b) < k, we must have 
n=l, i>(c) = k — 2 and t>(&) = — 1. Hence (3 = 1 mod 2. Moreover, 
if n = 1 we must have e f^fr) = e (|). Thus 



G(6,c)= £ e(^J=2* 

2 mod 2 fc ^ ' 

since + 2; = mod 2 for all z. □ 
Proposition 18. T7ie following equality holds for all x, y 
e (^(- fe2 + 2z (v ~ dx ) + ad 2 x 2 )J 



z mod 2 k 



z mod 2 k 
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Proof. The case v(b) > k, i.e. 6 = mod 2 k , implies that d = 1 
mod 2 fc and equality holds trivially. We may thus assume that t> (b) < k. 

We begin by noting that since y — dx = d(d~ 1 y — x) = d(y — x + aby) 
we see that v(y — x) < v{b) implies that v(y — dx) < v(b); putting 
x 1 = d~ x x we see that the converse holds, and hence v (y — x) < v(b) if 
and only if v(y — dx) < v(b). 

First case, v(y— x) < v{b): Putting c = y—x, c = y—dx respectively 
and applying lemma [TJ] we see that both sides are zero except when 
v(c) = k — 2 and v(b) = k — 1. For the exceptional case we note that 
v(b) = k — 1 implies that d^ 1 = 1 + ah = 1 mod 2 fc , and the same 
holds for d. Moreover, v(c) = k — 2 means that x = y mod 2 k ~ 2 and 
since 4|a we see that 

LHS = 2 k e ^ad 2 x 2 ^j = 2 k e {^ady 2 ^j = RHS. 

Second case, v(y — x) > v(b): As remarked above this means 
that v (y — dx) > v(b). We may thus complete the squares inside the 
exponentials, and we get 

LHS= £ e ( T - (-b(z - y —^) 2 + iV ~ h dX? + ad*x 2 



mod 2 k 



and 



z mod 2 fe 

After changing variables and taking constants outside we get 



LHS = e[^^ y - dxf 



2 k V b 



+ ad 2 x 2 )) £ eg(-fe 2 )) 

^ ^ z mod 2 fc 



and 



2 fe V b 



' ' ~ , ] ofc 



2 mod 2 fc 

Now, d = 1 mod 16 means that is a square modulo 2 fc . Changing 
variables by z — > ydz in the second sum we see that the sums are 
equal, and we are left to prove that 



2 k \ b / I \2 k 
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This will follow from the equality 

{V ~ u dxf + ad 2 x 2 = ^ y - x)2 + ady 2 . 
b b 

Collecting terms it is equivalent to 

= ad(y 2 - dx 2 ) + b^^dy 2 + dx 2 - 2dxy - y 2 - d 2 x 2 + 2dxy) 

= ad{y 2 - dx 2 ) + ft" V(d - 1) + x 2 {d - d 2 )) 

= ad(y 2 - dx 2 ) + (d- l)b-\y 2 - dx 2 ), 
which follows from the identity 

ad+ (d-l)/b = d(a+ - — -^) 

l-(l + a6). 

ab. 

= d(a-—) = 0. 
b 



□ 
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Appendix B. Counting norm one elements 

Let e be the ramification index of a prime p in Ok, i-e. the largest 
integer such that P s \pOk, where P C Ok is any prime ideal dividing 
P&K- Since K is quadratic e G {1,2}. If e = 2 then p is said to be 
ramified. If e = 1 then p is called unramified, and one of two things 
can happen: either pOx = P, in which case p is said to be inert, or 
pOx = PP, in which case p is said to split. 

Now, fix a prime p with ramification index e, be it one or two. The 
norm map 

Af : K —> Z 
descends modulo p k , and gives a homomorphism 

{o K /P ek Y - (z//) x . 

We let 

(d^/p 6 *) 1 

be the kernel of this map, i.e., the group of norm one elements. For 
I < ek we let 

((l + p'ytl + F*)) 1 

be the norm one elements in the subgroup (1 + P l ) /(l + P efc ), these are 
precisely the norm one elements that reduce to one modulo P l . 

Lemma 19. We have 

{{p — l)p k ~ l if p is split, 

(p+l)p k ~ 1 if pis inert, 

2p k if p is ramified. 

Proof. Recall first from class field theory ||] that the index (in Z*) of 
the image of the units in the p-adic completion of Ok under the norm 
map equals the ramification index e. We will split the proof in three 
parts: 

The split case: If p splits in K then pOx = P\Pi where P\,Pi are 
prime ideals in Ok, and where P2 = Pi- The map x — > x gives an 
isomorphism between Ok/P± and Ok/ Pi- This, together with the 
Chinese Remainder Theorem gives 

K /p k O K ^ Ok/P? x Ok/ Pi ^ Ok/ Pi x K /Pi 

where x G D^/p fc Dx is mapped to (x,x) G Ok/ Pi x Ok/ Pi- Fur- 
thermore, Ok /Pi — Z/p fc Z, and therefore 

(B.l) K /p k O K ^ Z/p k Z x Z/p k Z. 

Under this isomorphism, Galois conjugation maps (x,y) G Z/p k Z x 
Z/p k Z to (y, x). Thus the natural embedding of Z/p k Z in Di^/p fc Dx — 
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Zjp k Z xZ/p k Z consists of elements of the form (x, x), and the image of 
(x, y) under the norm map is (xy, xy). Hence the norm one elements in 
&K/p k &K correspond to elements of the form (x,y) G Z/p k Z x Z/p k Z 
such that xy — 1, and the number of such elements is (p — 

The inert case: Here e = 1 and the local norm map is onto Z*; 
reducing modulo p we get an exact sequence 

1 - {0 K /P k f - (Dk/PT - (Z//) x - 1. 

Hence 

K^)'l = ^gf = (?+!)/-'■ 
The ramified case: Here the image of the norm map in Z* is of index 

2, and thus the image of the norm in (Z/p k ) x has cardinality <yP ~ 1 ]^ — . 
Consequently, 

|(o^p efc )V ol(DWpefc)X| 



(p — i)p k 1 

Now, 

\(D K /P ek r\ = \(Dk/P) x \ x \(1 + P)/(1 + P ek )\ = (p-l)p ek ~ l 
and since e = 2 we get 

lP ek \ l \ = 9 

(p — l)p fc_1 



□ 



We will also need to know the number of norm one elements that 
reduce to one modulo P l . 



((i+p i )/(i+p ek )Y 



Lemma 20. We have 

p k ~ l if p is split or inert, 

K p x p k +\ l / 2 ~\- 1 if p i s ramified. 

where K p = 1 if p is odd, and K 2 = 1 or 2. 

Proof. The split case: From the previous discussion of the isomorphism 
in equation [B.l| we see that norm one elements congruent to one modulo 
P[ correspond to elements (x, x -1 ) G Z/p k Z x Z/p k Z such that x = 1 
mod p l . The number of such elements is |(1 +p l )/(l + p k )\ = p k ~ l . 

The inert case: If p is odd then x — > x 2 is an automorphism of 
(1 + P l )/(l + P k ) since the order of the group is odd. Thus the norm 
is locally onto in the sense that the map 

M : (1 + P l )/(l + P k ) -»• (1 + p l )/(l + p k ) 

is onto. 
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If p is even (and inert) then squaring is not an automorphism as 
(1 + x) 2 = 1 + 2x + x 2 . However, 1 + p l C 1 + P l and squaring maps 
(1 + p l )/(l + p k ) onto (1 + p l+1 )/(l + p k ). Thus 

(1 + p l+1 )/(l + p k ) c M ((1 + P l )/(l + P fe )) , 

which shows that the image of the norms must be either (1 +p l+1 )/ (1 + 
p k ) or (1 +p l )/(l +p k ). (There are no subgroups in between!) We will 
show that the former holds; since 2 is unramified the discriminant of 
K is odd and O k = Z[±±yj£]. Hence tr(0 K ) = Z, and there exists 
x E D K with odd trace. Now, 

M{\ + p k x) = l+p k tr(x) + p 2k M(x), 

shows that the image must be (1 +p l )/(l +p k ). 

Thus, whether p is even or odd, the norm map is locally onto and 
hence 

((l + ^/fl + P* 1 



((1 + P')/(1 + ^)) 1 



„k-l 



P 



|((l+pl)/(l +p fc))> 

The ramified case: First we note that 

(B.2) M((l + P l )/(1 + P ek )) C (l+p^)/(l+p k ). 

Arguing as before that squares are in the image of the norm we see 
that equality holds for odd p, and we obtain 

((1 + P')/(1 + P efc )) : 



((1 + /*)/(! + /**)) 



\0 K /P 



2k-l 



p 



((l +p R/2l)/(i +p fc)) ; 

2k-l 



pk-ll/2] pk-ll/2] 



p 



k+\l/2]-l 



For p even the squaring argument shows that 

(i + p ri/2i+i)/(i + c m ((i + p 1 )/(i + P ek )) , 

which gives a lower bound on the image. This gives the same result as 
for the odd case, except for a factor of 2. □ 



36 



PAR KURLBERG AND ZEEV RUDNICK 



References 

1. R. Balian and C. Itzykson Observations sur la mecanique finie, C.R. Acad. Sc. 
Paris 303 (1986) 773-778. 

2. F.R. Beyl The Schur Multiplicator of SL{2, Z/rnZ) and the Congruence Subgroup 
Property, Math. Zeit. 191 (1986), 23-42. 

3. A. Bouzouina and S. De Bievre Equipartition of the eigenf unctions of quantized 
ergodic maps on the torus, Comm. Math Phys. 178 (1996), 83-105. 

4. J. W. S. Cassels and A. Frohlich Algebraic Number Theory, instructional con- 
ference on algebraic number theory (1965 : University of Sussex, Brighton), 
Thompson 1967. 

5. Y. Colin de Verdiere, Ergodicite et fonctions propres du laplacien, Comm. Math. 
Phys. 102 (1985), 497-502. 

6. M. Degli Esposti Quantization of the orientation preserving automorphisms of 
the torus, Ann. Inst. Poincare 58 (1993), 323-341. 

7. M. Degli Esposti, S. Graffi and S. Isola Classical limit of the quantized hyperbolic 
toral automorphisms, Comm. Math Phys. 167 (1995), 471-507. 

8. S. Gelbart Weil's representation and the spectrum of the metaplectic group, Lec- 
ture Notes in Math. 530, Springer- Verlag, 1976. 

9. J.H. Hannay and M.V. Berry Quantization of linear maps on a torus - Fresnel 
diffraction by a periodic grating, Physica D 1 (1980), 267-291. 

10. D. Jakobson, Quantum unique ergodicity for Eisenstein series on 
PSL 2 (Z)\PSL 2 (R), Ann. Inst. Fourier (Grenoble) 44 (1994), 1477-1504. 

11. J.P.Keating, The cat maps: quantum mechanics and classical motion, Nonlin- 
earity 4 (1991), 309-341. 

12. S. Klimek, A. Lesniewski, N. Maitra and R. Rubin Ergodic properties of quan- 
tized toral automorphisms, J. Math. Phys. 38 (1997), 67-83. 

13. S. Knabe On the quantisation of Arnold's cat, J. Phys. A: Math. Gen. 23 (1990), 
2013-2025. 

14. T. Kubota On automorphic functions and the reciprocity law in a number field, 
Lectures in Mathematics, Department of Mathematics, Kyoto University, No. 
2 Kinokuniya Book-Store Co., Ltd., Tokyo 1969. 

15. W. Luo and P. Sarnak, Quantum ergodicity of eigenf unctions on PSL 2 (Z)\H. 2 , 
Inst. Hautcs Etudes Sci. Publ. Math. 81 (1995), 207-237. 

16. J. Marklof and Z. Rudnick Quantum unique ergodicity for parabolic maps, 



preprint math-ph/9901001 , to appear in Geom. and Funct. Analysis. 

17. J. Mennickc On Ihara's modular group, Invent. Math. 4 (1967), 202-228. 

18. A. Nobs Die irreduziblen Darstellungen der Gruppen SL 2 (Z p ), insbesondere 
SL 2 {Z 2 ), I. Comment. Math. Helvetici 51 (1976), 465-489. 

19. I. Percival and F. Vivaldi Arithmetical properties of strongly chaotic motions, 
Physica D 25 (1987), 105-130. 

20. Z. Rudnick and P. Sarnak The behaviour of eigenstates of arithmetic hyperbolic 
manifolds, Comm. Math Phys. 161 (1994), 195-213. 

21. A. Schnirelman Ergodic properties of eigenfunctions, Usp. Math. Nauk 29 
(1974), 181-182. 

22. I. Schur Untersuchungen iiber die Darstellung der endlichen Gruppen durch 
gebrochene lineare Substitutionen, J. Reine Angew. Math. 132 (1907), 85-137. 



HECKE THEORY AND EQUIDISTRIBUTION FOR QUANTIZED CAT MAPS 37 

23. O. Taussky Introduction into connections between algebraic number theory and 
integral matrices, appendix to H. Cohn A classical invitation to algebraic num- 
bers and class fields, Springer, New York 1978. 

24. S. Zelditch, Uniform distribution of eigenf unctions on compact hyperbolic sur- 
faces, Duke Math. J. 55 (1987), 919-941. 

25. S. Zelditch Index and dynamics of quantized contact transformations, Ann. Inst. 
Fourier (Grenoble) 47 (1997), 305-363. 

Raymond and Beverly Sackler School of Mathematical Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel. Current address: De- 
partment of Mathematics, University of Georgia, Athens, GA 30602, 
U.S.A. (kurlberg@math.uga.edu) 

Raymond and Beverly Sackler School of Mathematical Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel (rudnick@math.tau.ac.il) 



